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Abstract 

We studied the mathematical relations between species abundance distributions (SADs) 
and species-area relationships (SARs) and found that a power-law SAR can be gener- 
ally derived from a power-law SAD without a special assumption such as the "canonical 
hypothesis" . In the present analysis, an SAR-exponent is obtained as a function of an 
SAD-exponent for a finite number of species. We also studied the inverse problem, from 
SARs to SADs, and found that a power-SAD can be derived from a power-SAR under 
the condition that the functional form of the corresponding SAD is invariant for changes 
in the number of species. We also discuss general relationships among lognormal SADs, 
the broken-stick model (exponential SADs), linear SARs and logarithmic SARs. These 
results suggest the existence of a common mechanism for SADs and SARs, which could 
prove a useful tool for theoretical and experimental studies on biodiversity and species 
coexistence. 
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1. Introduction 



A fundamental question in ecology is how various species coexist in nature (Chave et 
al. 2002; Hutchinson 1959; Levins 1970; May 1972; Pacala k Tilman 1993; Rosenzweig 
1995; Tokeshi 1999; Gaston & Blackburn 2000; etc.). The answers to this question are 
expected to provide great insights into both theories of biodiversity and effective nature 
conservation practices. 

Among the various explorations into the mechanisms of species coexistence, two 
community-level properties have been theoretically and quantitatively examined: species 
abundance distributions (SADs) (Motomura 1932; Fisher et al. 1943; Preston 1948; 
MacArthur 1957; May 1975; Sugihara 1980; Harte et al. 1999; Hubbell 2001) and species- 
area relationships (SARs) (Arrhenius 1921; Preston 1962a,1962b; MacArthur & Wilson 
1967; May 1975; Pueyo 2006). A mechanism to reproduce such macroscopic ecological 
patterns using a microscopic model has been one of the central issues in recent community 
ecology (Durrett & Levin, 1996; Ney-Nifle & Mangel, 1999; Tokita & Yasutomi, 1999; 
Bastolla et al, 2001; Tokita & Yasutomi, 2003; Tokita, 2004; Lawson & Jensen, 2006; 
Tokita, 2006). Studies on these macroscopic patterns, therefore, not only give theoretical 
insight into large-scale ecosystems but also clarify the impacts of habitat fragmenta- 
tion. Thus, they aid in efforts to devise long-term estimations and strategies for nature 
conservation. 

SAD and SAR are mutually connected to each other. Preston (1962) derived the 
power-law SAR from the lognormal SAD under an assumption called the canonical hy- 
pothesis, which states that the peak of the individuals curve coincides with the num- 
ber of individuals in the most abundant species. May (1975) comprehensively studied 
various types of SAD such as Preston's lognormal distribution, uniform distribution, 
MacArthur's Broken Stick (exponential) distribution, Motomura's geometric series dis- 
tribution and Fisher's logseries distribution, and demonstrated that the first three SADs 
lead to power-law SAR and the latter two correspond to log SAR. 

In addition to those pioneering works on SADs, power-law SAD has been reported 
(Margalef 1994; Pueyo 2006). Other than SADs, power-law has been observed for rela- 
tionships between abundance and body size (Siemann et al. 1996). Power-law is, there- 
fore, ubiquitous in biology. It is moreover known that, in general, power-law distribution 
is given in a limit of large variance of the lognormal distribution. 

Here, we demonstrate that power-law SAR can be mathematically derived from 
power-law SAD without any such assumption as the canonical hypothesis. We also 
discuss an inverse problem: namely, what type of SAD can be obtained when we start 
from the power-law SAR? 

2. Species abundance distribution and rank 

Let a(x)dx be the species abundance distribution (SAD) between the number of individ- 
uals x ~ x + dx. The total number of species S is obtained by integrating cr(x) from the 
minimum value of x = m to the maximum value of x = X as 




(1) 
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The species rank of the number of individuals x is defined by 

POO 

R(x) = / a{x')dx'. (2) 



The inverse function of eqn ©, xr, is the rank-abundance distribution; that is, the 
number of individuals of the -R-th rank is xr. The first rank of species, i.e. the most 
abundant species, has the maximum number of individuals, X; then we obtain 

R(X) = 1. (3) 

This equation (j3J) is equivalent to the estimation used by Preston (1962a, 1962b), and 
May (1975). 

Using the SAD a(x), the total population or the total number of individuals becomes 

f x 

N = xa(x)dx. (4) 



Writing the number of individuals normalized by the minimum value m as x = x /m, N = 
N/m, and a(x)dx = a(x)dx, we obtain 

rX rOO pX 

S = I a(x)dx, R = / a(x')dx', N = xa(x)dx. 

Jl Jx Jl 

Hereafter, we omit the hat and we consider x, X, a, and N to be the normalized 
quantities. 

3. Power-law SAD 
3.1. Power-law SAR 

We consider here an SAD which decays with power-law from the minimum number of 
individuals, x = 1 (the number normalized by the minimum value), to the maximum 
value, x = X as 

a{x) = Sax-^ 1+a \ (5) 

where S is a constant. The rank function R{x) is calculated by integrating u{x) from x 
to infinity 

R(x) = / a(x')dx' = Sa = Sx^. (6) 

Jx -a x 

Using this, we find 

^oo rX poo 

S= a(x)dx = / a(x)dx + / a(x)dx = S + 1. (7) 
Ji Ji Jx 

where eqns © and © are used for the last equality. 

The relation of the total number of species S and the maximum number of individuals 
X is obtained by eqns © and ©, 1 = R(X) = SX~ a , so we obtain 

X a = S = S + 1. (8) 
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By inserting eqn (JSJ) into eqn @, the total number of individuals becomes 

rX j£ X ' 

N = xa(x)dx = a . (9) 

Ji 1 - a 

For very large X, from eqn (jHJ), N ~ aX/(l — a) and S oc N a for a < 1, and iV ~ 
aX a /(a — 1) and S oc N for a > 1. When the size of area A is proportional to the 
normalized total population N, the SAR becomes S = cA z with z = a for a < 1 and 
2 = 1 for a > 1. 

For a finite value of X, the SAR never becomes a simple power daw relation. In this 
case, let us start to define the SAR-exponent £ as 

dlnS" . . 

c =dbX (10) 

This can quantify the increasing rate of the species richness with increasing area size A. 
For the power SAR, S = cA z , this exponent becomes £ = z, and for the logarithmic case, 
S = K Yn.A + a, £ = l/(ln A + a/K). This exponent is closely related to the persistence 
function a(A) introduced by Plotkin et al. (2000) as £ = — loga(^4). 

If the normalized total population ./V is proportional to the area size A, this SAR 
exponent equals £ = dlnS/dln 2V. In the case of the powerdaw SAD eqn ©, we obtain 

C = z!l = — = (ii) 

<> d/U 1 + 0_ e -fly g ' 

+ i _ e -e y 



where /i = In 2V, y = In S ~ In S using S ~ <§ S> 1 for a large community, and = — — 1. 
For large N and A and in the limit of an infinite number of species (S — > oo), ^ becomes 
constant: 

I a (for a < 1) . . 

C = < V 7 . (12) 

[1 (faro>l) 

In Fig. 1, we show the SAR exponent £ v.s. a. The £ bends near a = 1. 



Figure 1: SAR exponent £ v.s. the exponent a for the powerdaw SAD eqn (jHJ) 
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3.2. Logarithmic SAR 



For a — > 0, keeping Sq = Sa constant, the SAR becomes logarithmic, because for the 
SAD, a{x) = S ar( 1+a ) and 

x 



S 



N 



x 



a(x)dx 



So In X, 

x 

xa{x)dx 



S x~ a 



a 



So 



l-X~ 



Sr 



d a (l — exp(— a lnX)) 



o 



«=o 



So 



x 



x a dx 



- a 



-{X 



l-a 



S X (for a -►()), 



(13) 



(14) 



where d a denotes the derivative in a, and L'hopital's rule is used in taking the limit in 
eqn ED then we derive 

f~mf (i5) 

bo Do 

This case a — > corresponds to the continuous version of the geometric SAD (Motomura 
1932, May 1975). In the geometric SAD, the rank-size distribution is X{ = NCkk^l — ky -1 
with a constant k, and the coefficient Ck being the normalization constant given by the 
condition J^^ =1 x\ = N, = 1/[1 — (1 — k) s ]. The inverse function of this expression of 
Xi leads to the rank function by setting i = R, 

i = R(x) = l + K ln(NC k k) -Klnx 

Differentiating R{x) by x, we obtain 

a(x) 



K 



In 



dR(x) 
dx 



K 

x 



(16) 



(17) 



If we express it by the normalized quantities, X = x\ = NCkk, 1 = xs = NCkk(l — k) s . 
For large N, S, and X, the parameter Ck becomes Ck — ► 1, and X ~ kN, S = K\n(kN); 
finally, we obtain the logarithmic SAR. From eqn (|17[). we find that the SAD corresponds 
to the power-law SAD eqn Q with a — * 0, and we can obtain the logarithmic SAR for 
the power-law SAD with a — > in the continuous approximation as well. The relation 
between parameters Sq, K and k is So = K ~ 1/k. 



4. From SAR to SAD 

In contrast to the previous section, we consider here the inverse problem of what kind 
of SAD can be derived from a given SAR. If the shape of the SAD is unchanged with 
increasing S, X, and N, we can write cr(x) = Sp(x): only the coefficient S varies with 
varying S, but the function p(x) is unchanged. In this case, eqn (^Q) becomes 

S = sj^ p(x)dx, (18) 

with the use of the normalized quantities (normalized by the minimum number of indi- 
viduals). Defining the cumulative distribution function as 

P{x) = J™ p(x')dx' (=^f), (19) 
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we obtain 

1 



s , P(X) = I p(x)dx, (20) 



x 



from eqns © and ©. 

From eqn @, the normalized total population N divided by S is 



N rX 



f 

/ xp{x)dx. (21) 
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If p(x)dx = 1, from eqns (|T8|) and (|20j) . we obtain 

roo pX foo 

S = S p(x)dx = S p(x)dx + S p(x)dx = S + 1 ~ S. (22) 
Ji Ji Jx 

Equations (|2()|) is the relation between X and S through P{X) on the one hand, and eqn 
(|21|) is the relation between and 5" on the other. 

The SAR is the relation between N and 5 if oc A, and if we put A^ = SF(1/S), 
from eqn (J20[). A^/S = F(P(X)). Therefore, if X is changed to X + OT, the variation of 
eqn (|2Tj) becomes 



5F(P(X)) = -X5P(X) + I x dp ^ dx6/3, (23) 

where j3 expresses some parameters of the distribution, and we use SP(X) = —p(X)5X. 
Because we consider the case in which the parameters are unchanged with varying N, S, 
and X, the last term vanishes. For example, if the SAR is power-law S = cN z , N = 
c 'gi/z _ (J p(x)~ l l z ; that is, F{P) = c'P l ~ l l z . Substituting this expression into eqn 
(|2*3|) . we obtain 

c' (l - -^j P(X)- 1 /^P(X) = -X5P(X). (24) 

Therefore, the exponent z must satisfy 

z < 1 and P(X) oc AT 2 . (25) 

This means that the tail of the SAD is a power-law of eqn and the exponent becomes 
a = z < 1. 



5. Discussion 

We obtained the power-law SAR for a power-law SAD using the classical method of R. 
May and Preston. We also considered the inverse problem of obtaining the power-law 
SAD from a given power-law SAR if the shape of the SAD is unchanged with varying 
total population, area size, and species richness. 

R. May (1975) obtained the power-law SAR for the lognormal SAD and the broken 
stick model SAD. In both cases, the parameters of SAD vary with increasing A^: in the 
lognormal SAD, the variance of the log of the number of individuals increases proportional 
to In AT, and for the case of the broken stick model, in which the SAD becomes an 
exponential distribution, the SAR is obtained if the average of the number of individuals 
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is proportional to the total number of species S. The last term of eqn Q23|) does not 
vanish in either SAD. 

In the case of a linear SAR, S = cN, first we consider that this linearity holds in 
the limit of S — ► oo and assume that any deviation from the linearity becomes N/S = 
1/c - bS~^(j > 0), then from eqn J23J) IrfP^SP = X5P and P(X) oc X' 1 ^ 1 '^. 
Hence, this case corresponds to the power-law SAD with a = 1/(1 — 7) > 1. If the 
linearity holds completely for all S as b = above, the variation of the maximum number 
of individuals becomes zero, 5X = 0, and this corresponds to the case in which X is 
constant, for example, a uniform SAD (May 1975), cr(x) cx 5(x). 

We found in the previous section that the logarithmic SAR, S = K In N+a, is derived 
from the case of the SADoc 1/x, and the rank function defined by eqn Q diverges, so 
we cannot use eqn (|23|) as a method for deriving SAD from SAR. 

Harte et ai. (1999) obtained an SAD from a power-law SAR using the renormalization 
group technique for the case of existing self-similarity: the fraction of a species found in 
an area with a size A, which is also found in A/2, is independent of A and the abundance 
x. But this SAD is far from the power-law SAD. Pueyo (2006) pointed out that there 
are possibilities for other shapes of SADs if the fraction depends on the abundance, and 
he also obtained power-law SAD from more straightforward discussion. 
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